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ABSTRACT 

Observations of clusters in the 30 to 350 GHz range can in principle be used to determine a galaxy cluster's Comp- 
tonization parameter, y, peculiar velocity, v and gas temperature, Tg via the dependence of the kinetic and thermal 
Sunyaev-Zeldovich (SZ) effects on these parameters. Despite the significant contamination expected from thermal 
emission by dust in high-redshift galaxies, we find that the simultaneous determination of r, v and Tg is possible 
from observations with sensitivity of a few fiK in three or more bands with arc minute resolution. After allowing for 
realistic levels of contamination by dusty galaxies and primary CMB anisotropy, we find that simultaneous determi- 
nations of velocities to an accuracy of better than 200 kms^^ and temperatures to roughly keV accuracy should be 
possible in the near future. We study how errors change as a function of cluster properties (angular core radius and 
gas temperature) and experimental parameters (observing time, angular resolution and observing frequencies). Con- 
taminating synchrotron emission from cluster galaxies will probably not be a major contaminant of peculiar velocity 
measurements . 

Subject headings: 



1. INTRODUCTION 

The cosmic microwave background (CMB) is a tremen- 
dous tool for studying cosmology. The anisotropics in the 
CMB provide a wealth of cosmological information (for 
a recent review see Hu & Dodelson 2002) and upcoming 
experiments will provide high sensitivity (approaching 1 
fj.K) and high angular resolution (approaching 1'). Such 
precise measurements will provide interesting constraints 
on secondary anisotropics in the CMB, imprinted by ma- 
terial along the line of sight at redshifts z ^ 1000. 

The largest secondary anisotropy is expected to 
be caused by galaxy clusters, through the S u nyaey - 
Zeldovich (SZ) effect (ISunvaev fc Zel'dovichl IiQTI . 
Compton scattering of CMB photons by intra-cluster 
electrons leaves a spatial and spectral imprint in the 
CMB. Recent measurements of the SZ effect have pro- 
vided detailed maps of the electron distribution in galaxy 
clusters (see Birkinshaw 1999 and Carlstrom et al. 2002 
for recent reviews) but new experiments offer the promise 
of an order of magnitude increase in sensitivity. With 
sensitivity it should be possible to probe other cluster 
properties, such as bulk velocities and electron tempera- 
tures. 

In this paper we investigate the feasibility of measur- 
ing peculiar velocities and cluster temperatures at cm 
and mm wavelengths in the presence of contaminating 
sources, focusing on the effects of dusty high-redshift 
galaxies. These iss ues have been partly addressed pre- 
viously bv .Fischer fc Langd l(T993.) and .Blain UMSj but 
both our understanding of dusty distant galaxies and 
CMB experimental capabilities have advanced tremen- 
dously in recent years. 

We forecast errors by calculating the SZ-parameter 



Fisher matrix for multi-frequency maps that contain 
signals from the clusters as well as noise from CMB 
anisotropy, dusty galaxies and the measurement pro- 
cess. Our work exten ds previous studies of SZ 
parameter reconstru ction ( Aghanim "eTaOllOOll l200a 
IHaehnelt fc Tegrnark.l99& : .Holder ,2003i) by considering 
point source and CMB anisotropy contamination and the 
influence of the atmosphere on the relative sensitivity of 
different frequency bands. While Haehnclt & Tcgmar^ 
(^,996) did include the effect of CMB anisotropy con- 
tamination, explicitly subtracting it with an optimal fil- 
ter, point sourc e conta mination was ignored. And while 
lAghanim et all l|200H) included both CMB and point 
source contamination, forecasts were for one specific ex- 
periment: Planck. Here we vary experimental parame- 
ters in order to guide experimental design. 

In §2 we outline the SZ effect and how it relates to 
cluster properties. §3 summarizes the relevant known 
properties of dusty galaxies observed in the direction of 
galaxy clusters, while §4 outlines our methods for esti- 
mating the effects of contaminating sources on cluster 
measurements. We discuss complications that will arise 
in analyses of real data in §5 and describe the experi- 
ments we consider in §6. Our primary results are in §7, 
where we forecast uncertainties for a four-channel refer- 
ence experiment as a function of cluster gas temperature 
and core radius. We explore the impact of varying exper- 
imental parameters such as angular resolution, observing 
time and number and placement of frequency channels. 
We follow that with a discussion of possible radio point 
source contamination and close with a discussion of our 
results and implications for future instrumentation. 

2. SUNYAEV-ZELDOVICH EFFECTS 
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The SZ effect is the change in energy of CMB photons 
from Compton scattering with electrons, pri marily in the 
intra-c luster mediu m. For recent re views see iBirkinshawl 
ljl999H and Carlstr om et all lj2002() . The high tempera- 
ture (several keV) of the electrons relative to the CMB 
photons leads to a net increase of energy of the scat- 
tered photons, while a bulk motion of the cluster elec- 
trons along the line of sight leads to either a net red- 
shift (moving away from observer) or blueshift (moving 
toward observer) of the scattered CMB photons. The 
first effect is the thermal SZ effect and it has a distinctly 
non-thermal spectrum; the scattering preserves photon 
number and mainly just shifts the energy of each pho- 
ton, distorting the original blackbody spectrum. The 
effect of the bulk motion is the kinetic SZ effect, and the 
emergent spectrum is that of a blackbody with a slightly 
different temperature. 

For the thermal SZ effect, the main physical parame- 
ters are the optical depth to Thomson scattering, r and 
the fractional energy gain per scattering O = kT^ / m^.c^ . 
The combination tQ sets the amplitude of the spectral 
distortion. For the kinetic SZ effect, the redshift or 
blueshift is set by /3 = ?;/c, so the relevant combination 
is t/3. 

Typical electron temperatures are on the order 
of 0.01 me^\ making relativistic effects modestly 
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) . The relativistic corrections are sen- 



sitive to the electron temperature through the relativistic 
correction to the Thomson cross-section and through the 
relativistic form of the thermal Maxwellian velocity dis- 
tribution. These corrections can be on the order of 10% 
or higher at many observing frequencies. 

Thus the temperature difference due to the kinetic SZ 
effect is given by 

^=rP (1) 

cnib 

while the temperature difference due to the thermal SZ 
effect as a function of frequency v (expressed in dimen- 
sionless units x = hv/kTcmb) is given by 



cmb 



-4)Te[l + ^(x,Te)] 



(2) 



where the relativistic corrections have been represented 
by d{x,Te). 

Typical values for massive clusters ofr ~ 0.01, 6 ~ 
0.01 and P ~ 0.001 lead to typical amplitudes (AT/Tcmb) 
for the thermal effect, kinetic effect and relativistic cor- 
rections of 10"^, 10~^ and 10~^ respectively. 

In this work, we take our fiducial galaxy cluster to 
have V = -200km s-\ = 6 keV and r = 0.01. We 
vary th e thermal SZ propertie s using the scaling rela- 
tions of iMcCarthv et all l|2003H . The central y parame- 
ter was taken to scale as y oc T^, which leads to r oc T. 
Self-similar evolution would result in y cx T'^/^; we use 
the slightly steeper relation to be consistent with evi- 
dence for excess entropy. We ignore evolution of intrin- 
sic cluster properties with redshift, such as the expecta- 
tion that clusters (of a fixed mass) become more com- 
pact (and therefore hotter) at higher redshift. It is ex- 
pected that the effect of this evolution is to enhance the 



thermal SZ effect, but the details of this redshift evo- 
lution are sensitive to poo rly understood gas processes 
lIHolder & Carlstromll200l . 

3. MILLIMETER EMISSION FROM DUSTY GALAXIES 

Star-forming galaxies can be very luminous in the sub- 
millimeter wavelength range, as recent measurements 
by the Sub-millimeter Common User Bolometer Ar- 
ray (SCUBA) and other experiments have discovered 
(see Blain 2002 for a recent review). In particular, 
galaxy clusters are often found to coincide with relatively 
bright submillimeter er nitters, probablv b ecause of gravi- 
tational lensing effects l|Smail et al.ll200^ . The emission 
is mainly from dust at a temperature of several tens of 
Kelvin, putting the peak of the radiation spectrum at 
submillimeter wavelengths, but there can still be signifi- 
cant emission at millimeter wavelengths. 

Many galaxy clusters have been observed to coincide 
with submillimeter sources with fluxes of several mJy 
at ly—SbO GHz. Assuming a typical spectrum for these 
sources would lead to expectations of fluxes at 150 GHz 
on the order of 1 mJy. In a beam of I'xl' this would 
correspond to nearly 30 fiK of contamination if not cor- 
rectly taken into account. This makes dusty galaxies a 
non-negligible contaminant for studies of galaxy clusters 
at mm wavelengths (Blain 1998; Fischer & Langc 1993). 

Throughout we often characterize a given intensity as 
the equivalent departure from the mean CMB temper- 
ature, AT. The conversion factor, for AT << Tcmb, 
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3.1. Shot Noise 
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The experiments we consider are not sensitive to the 
absolute flux, but to spatial variations in the flux. These 
variations contribute to the variance of temperature fluc- 
tuations on the sky, observed with a Gaussian beam pro- 
file with full width at half maximum of 0h, 
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f is the 



This variance decreases with beam size as AT^ cx 1/0^. 
For the convenience of using a quantity that is a property 
of the sky only, and not the angular resolution of the 
telescope, we define 5T^ = AT^ilf, where fti, oc 
solid angle of the beam. 

The shot noise at 350 GHz can be inferred from 
SCUBA observations. Borvs et al. (20031) used these 
observations to fit a double power-law form for dN(> 
S)/dS: 

s-A" f S ' ' 
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with = 1.8 mJy, Nq = 1.5 x lO^deg^^ a = 1.0 
and b — 3.3. The result is shot noise with an rms of 
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= 170/iK-arcmin where the subscript "350" is used 
to denote the frequency. 

3.2. Spectral Variations 

If the spectral behavior of these sources is known, then 
it would be straightforward to measure the flux at a 
higher frequency, where the emission is stronger and the 
SZ effects are smaller, and subtract the appropriate levels 
from the measurements at other frequencies. However, 
the spectral behavior is not perfectly known, so there will 
be uncertainty due to imperfect subtraction, as well as 
uncertainty due to noise in the measurement at higher 
frequencies. 

To investigate the homogeneity of the spectral behav- 
ior, we used submillim eter observations of local galax- 
ies l)Dunne et al.l l2?)(¥it) . Local galaxies were modeled 
as having a spectrum of the form / cx v^^'^B^, where 
BuiTdust) is the blackbody function and Prnm here refers 
to the emissivity spectral index (not the bulk velocity of 
the cluster!). Dunne et al. provide best-fit values of Pmm 
and dust temperatures. We used these best-fit spectra 
for the local galaxies, placed them at a range of redshifts 
and fit them over the range 150 to 350 GHz assuming a 
power-law (x v". As a function of redshift the mean 
spectral index was well-fit by a = 3.14 — 0.22z, with 
a fairly constant scatter at each redshift of ctq, ~ 0.18. 
The moderate flattening with increasing redshift arises 
because the radiation that is observed at 350 GHz origi- 
nates at a higher frequency, closer to the peak of the dust 
emission. 

Without source redshifts, it is not clear what to take as 
the mean spectral index. Assuming a uniform distribu- 
tion in redshift between and 5, the mean spectral index 
is 2.6 and the rms scatter around this mean is 0.4. The 
sources are unlikely to be at very low redshift, simply due 
to volume considerations, but the redshift distribution of 
the observed submillimeter-emitting galaxies is very un- 
certain. We set our mean value for a as a = 2.6 and set 
ctq = 0.3, unless specified otherwise. We set aa < 0.4 
since the redshift distribution will be more peaked than 
the uniform one that leads to 0.4. 

The quantity aa is the spectral index rms for individ- 
ual galaxies. For N galaxies of similar brightness the 
composite intensity will have a scatter in spectral index 
that is smaller by ^/N. We define an effective number 
density of galaxies by weighting them according to their 
contribution to shot-noise variance. Thus, 
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dlnS 



dS 

= 0.64/(sq.arcmin) (7) 

where the last equality follows from the iBorvs et al.l 
l|2003(l dN{> S)/dS. For simphcity we set N^s = l/(sq. 
arcmin). Finally we define (Tq, = Ca/VNcs and take our 
fiducial value to be CTq, — 0.3-arcmin. 

3.3. Gravitational Lensing 

dN{> S)/dS is for an unlensed 



The lBorvs et al.l l|2r 

population of sources. Lensing changes both the bright- 
nesses of galaxies and their number densities. Indeed, 
the magnification provided by galaxy clusters has been 
exploite d to study the population to fainter limiting mag- 
nitudes ijSmail et al.lll997l I200I) . The magnitude, and 



even the sign, of the effect on shot noise depends on 
dN{> S)/dS. Magnification increases S which increases 
shot noise, but also decreases dN{> S)/dS (at fixed un- 
lensed S) which decreases shot noise. For a power-law, 
dN{> S)/dS (X S"', the net effect is an increase in shot 
noise if a < —2. For the relatively bright population at 
350 GHz of interest for this work, the shape of dN/ dS is 
such that lensing leads to an enhancement of shot noise. 

Blain (1998) pointed out that lensing would signifi- 
cantly enhance the level of contamination by IR luminous 
galaxies of SZ observations. The key quantity for lensing 
is the Einstein radius, 0e- A point source with source 
plane location coincident with the cluster center {9 = 0) 
will form a ring of infinite magnification appearing in the 
image plane at distance 9e- In the image plane, magni- 
fication rises from the center towa rds 9e and then drops 
beyond 9e- Thus, as l)BlaiiJll9'9^ showed, the enhance- 
ment of confusion noise for an observation towards the 
cluster center increases with increasing beam size until 
peaking at 9^ — 29e as more weight is placed on galaxies 
near 9e, and then drops with further increase in 9i, as the 
effect is diluted by adding in more galaxies from beyond 
the Einstein radius. 

For a rough model of the effect of lensing we introduce 
the confusion noise enhancement factor, E{9b,9) which 
is the ratio of lensed confusion noise to unlensed confu- 
sion noise at distance 9 from the cluster center and for 
observations with angular resolution 9h - For t he dep en- 
dence of E on 9b we use the results of iBlainI l)1998f) for 
his galaxy evolution model II for a cluster with velocity 
dispersion of — 1360 km/s and 9e — 16" when it is 
placed at z = 0.171. For this case E{9b,0) = 1.1, 1.5, 
2.5, 2.1, 1.8 and 1.5 for 9b = 10, 20, 40, 60, 100 and 200 
arcsec respectively. We ignore the redshift-dependence of 
E since for 9b > 30" lBlaii]| l|1998^ shows that E is highly 
independent of the cluster redshift. Wc also ignore the 
<7y dependence of E since iBlainI l|1998fl shows this de- 
pendence to be quite slow, covering a range of 30% as 
ay varies from 800 km/s to 2000 km/s. For ay — 1360 
km/sec one expects a gas temperature of about T = 10 
keV (Girardi ct al. 1996). For dependence on 9 we sim- 
ply set E{9b,9) = E{9b,0) for 9 < 9b and E{9b,9) = 1 
otherwise. For all cases that we consider here the beam 
size will be larger than a typical Einstein radius for a 
cluster. 

4. METHODS 

We model the data at sky location 9i and frequency I'a 
as having contributions from both signal and noise: 

AT 

Tf: {Va,9i) = s{l^a,Oi) + n{l^a,9t). (8) 

JCMB 

We consider each component in detail below. 

4.1. Signal Model 
The signal is due to three components, 
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) = P-mh{Xi)fm{ya), (9) 



where 
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Fig. 1. — Derivatives of the frequency— dependent temperature 
fluctuation data, AT, with respect to the parameters Pi = y = t9, 
Pa = T/3, Pg = r6l2, P4 = AT350/TCMB and P5 = a (see text). 
The P3 derivative has been divided by 5 to fit on the plot and the 
P5 derivative is arbitrarily normalized. 
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dlaTTigV / c = T(3r and 



(11) 
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Here rig is the electron number density and = 
kTe/ {rrieC^) is the dimensionless gas temperature. The 
subscripts on the average quantities indicate the weight- 
ing along the line of sight. For r the weighting is with 
the electron-number density and for y the weighting is 
with the pressure, 71^0. The optical depth, r is given by 
T = / dlorne- 

The integrals defining Pi , P2 and P3 are along the line 
of sight through the cluster center. We take the angular 
profile of the cluster to be 



h{e)^{i + i\e\/e,y 



-1/2 



(13) 



The beam-convolved cluster profile is denoted by h{9), 
and we assume Gaussian beams in all that follows. Note 
that the assumption of the same cluster profile for all 
three signal components in equation El is only valid for 
an isothermal intracluster medium. 
Th e frequency dependences are given by ijltoh et al.l 



where 



(14) 
(15) 
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i = xcoth(x/2) 
s = 2;/sinh(a;/2). 



(17) 
(18) 



The first and third components are due to the scattering 
of CMB photons off of the hot electrons in the cluster and 
represent the thermal SZ effect, while the second com- 
ponent is the kinetic SZ effect. We have neglected con- 
tributions to the SZ effects that are t times higher order 
products of /3 (~ 10^'^) and Q 10^^). For clusters with 



temperatures of 10 keV or lower the higher order cor- 
rections are less than five percent llChallinor fc Lasenbvl 
119981 lltoh et alJll99a IStebbinsllT997|) at all freauencies 
except near the null of the thermal SZ spectrum (where 
the thermal SZ signal is sub-dominant). 

Because they are the amplitudes of the different fre- 
quency shapes. Pi, Pi and P3 are the theoretical param- 
eters most directly related to the data. From them one 
can estimate the physically interesting quantities: 



e, = P3/Pi = e, 



92 



(19) 
(20) 

(21) 



Thus the temperature determinations from SZ measure- 
ments are measurements of the pressure-weighted aver- 
age gas temperature. Velocity measurements are of the 
mass-weighted average velocity (assuming electron den- 
sity traces mass) times a correction factor dependent on 
differences in different temperature averages. The optical 
depth determinations are of the optical depth times an- 
other ratio of different average temperatures. Through- 
out, we assume isothermality in order to avoid these com- 
plications; the distinctions must be kept in mind when 
analyzing real data, particularly for comparison of X-ray 

and SZ gas temperatures; 

In a recent preprint ijHanse ^ l200l it was noted that 
peculiar velocity measurements of non-isothermal clus- 
ters could be biased. The source of this bias can be 
clearly seen in the above expression for where the true 
P is multiplied by the ratio of two different weightings of 
the electron temperature. In general the optical depth 
weighted temperature will not be the same as the pres- 
sure weighted temperature, resulting in a biased velocity 
estimate. This is a direct consequence of the biased es- 
timate of the optical depth. 

4.2. Noise Model 

The noise has a contribution from measurement error, 
riinst, from primary CMB fluctuations, ticmb and from 
emission from galaxies in the direction of the cluster rigai. 
We assume the composite dusty galaxy spectrum has a 
power-law form in intensity, Ii, cx v" where a is spatially 
varying with a{9) = a -I- da{0) and a — 2.6 as discussed 
earlier. Linearizing the dependence of on 6a we can 
write ritrai as 
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^)[fM + Sa{0,)M'ya)]E{0„0b) 

(22) 



where E{0i, 0^) is the shot noise enhancement factor dis- 
cussed earlier. The frequency dependences are given by 



/5(i^) = ln 



(e^ - l)V/4(350 GHz) (23) 
/4(^) (24) 



.350 GHz. 

The CMB noise contribution is independent of fre- 
quency, and the instrument noise contributions as a func- 
tion of frequency depend on the experiment being con- 
sidered. The total noise in the map is then ntot — 
ncMB + ^gai + ^inst- The noise in the map therefore has 
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the two-point function, considering points in the map 9i 
at frequency Va and 9j at frequency i>i,: 

ci:;r^Mer,va)ntot{ej,u,)) 
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To calculate errors on a new set of parameters, P, that 
are functions of our original parameters P, we Taylor 
■ffxp&iid the dependence of P on P to first order about 
"me fiducial value. We then calculate the Fisher matrix 
in these new coordinates by: 



<( 



In the above, c^^ is a shot-noise covariance ma- 
trix (with off-diagonal correlations only due to beam- 
smoothing) normalized to unit variance, and the last 
term gives the contribution of instrument noise, which 
will be considered in detail in §6. C'~^^^ is the covari- 
ance matrix of the CMB fluctuations, given by 

C™B ^ ^c,PKcos%)exp [-P0g/(21n2)] 
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(26) 

where 9ij is the angular separation between pixels i and 
J- 

This analysis implicitly assumes that we know the sta- 
tistical properties of all sources of noise, including the 
dusty galaxies. This is an excellent approximation for 
the CMB fluctuations, but the dusty galaxy noise power 
spectrum will be affected by both clustering and gravita- 
tional lensing. On arcminute scales the clustering effects 
should be relatively small, but the effects of gravitational 
lensing will modify the noise properties. 

4.3. Error Forecasting 

We forecast how well the three parameters of our model 
can be measured by calculating a Fisher matrix. Com- 
bining the spatial and spectral indicators i and a into a 
combined index fi (and j and b into v) we can write the 
Fisher matrix as 



ere the transformation matrix Rip — dPp/dPi. For 
example, such a variable transformation is necessary to 
get errors on = P3/P1. Dropping the higher-order 
terms in the Taylor expansion makes this procedure ap- 
proximate as well. 

Calculating these Fisher matrices can be computation- 
ally expensive as matrices of size Npi^Ni, have to be in- 
verted and multiplied where A^pix is the number of map 
pixels and is the number of frequencies. We therefore 
wish to pixelize no more finely than necessary, and keep 
the physical area of the map as small as possible. If the 
map is too small compared to the beam-convolved cluster 
profile we will lose the large-scale information necessary 
for subtracting off the CMB contamination. If the pixel 
size is too large we will lose small-scale information. We 
have found the following prescriptions to ensure we are 
neither losing information, nor using much larger A'pix 
than necessary: 



F , — 

-Tpp' — 



1 dAT, 



-'CMB 



dPr> 



(C"° 



.-1 dAT, 
't'" dPy 



(27) 



p ' ^^p' 

To this Fisher matrix we usually add a prior Fisher ma- 
trix which carries the information we have that is not 
contained in our measurements in the 30 to 350 GHz 
range. The error covariance matrix for our three param- 
eters is then given by 



pp' 



(F + Fp 



pp' 



(28) 



Note that since our data depend linearly on our param- 
eters P, this Fisher matrix calculation of the expected 
error covariance matrix is exact, given our model of the 
data. 

We include prior information due to X-ray measure- 
ments of the gas temperature with error ctt^ — ag x 
511 keV so that 

-2 



Pprior(l, 1) = (o-eT)" 



prior 
-fprior (3, 3) = 

r(l,3)=PpHor(3,l): 



: (aey) ^ and 



(29) 



Ppiioi (l; l)Pprioi (3, 3). 



Including the prior information in this manner is approx- 
imate since it results from Taylor expanding the depen- 
dence of 9 on Pi and P3. 



M = H29, + 9b) 
>c + 9b)/2.5 



■'pix 



(31) 

where the map and pixels are squares of length 9m and 
9pix respectively. 

5. POSSIBLE REAL-WORLD COMPLICATIONS 

Our model is fairly detailed and sufRcient in many 
ways, but there are several respects in which real clus- 
ters could present some challenges, such as temperature 
structure in the intra-cluster medium, incomplete knowl- 
edge of the lensing properties of galaxy clusters, and in- 
ternal bulk flows in the cluster. 

Temperature gradients in the intra-cluster medium 
would lead to a mis-match between the signal maps of 
equation Observed c lusters show sign ificant depar- 
tures from isothermality ijDe Crandi & Mo lcndi 2002.) at 
relatively large radii and at small radii, but the largest 
scales are already obscured by the primary CMB anisot- 
ropics and the small scales contribute relatively little to 
the SZ effect. It is therefore not expected to be a signif- 
icant effect but will complicate data analysis. 

Similarly, we have assumed that we know the cluster 
template. In the cases where we assume complementary 
X-ray spectroscopy it is quite reasonable to assume that 
the X-ray image gives an estimate of the electron spatial 
distribution. Without X-ray information the best esti- 
mate of the spatial template will come from the thermal 
SZ map itself, where a large mismatch between the as- 
sumed spatial template and the observed emission should 
be evident. A parametrized model for the cluster could 
be introduced into the fit without significantly affecting 
the constraints on peculiar velocities. 

We have assumed that the magnification due to lensing 
is a simple step function in radius. The main effect of 
lensing is to amplify the Poisson noise, and the amount of 
amplification depends on the details of the mass profile. 
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Incomplete knowledge of the mass distribution will there- 
fore lead to less efficient component separation. Deep 
optical images will be useful for the purpose of studying 
the strong lensing properties of galaxy clusters. 

On a related note, we have assumed that we know the 
statistical properties of the galaxy contamination. Lens- 
ing modifies the noise properties of the background galax- 
ies, and there is uncertainty in the statistics of dusty 
galaxies. Multi-frequency observations of many fields, 
both with and without clusters, will lead to a good un- 
derstanding of the statistics of dusty galaxies. 

We have also assumed a single peculiar velocity for 
the cluster whereas real clusters show evidence for in- 
ternal flows, sometimes as large as 3-4000 k m s~^ 
(jDupke fc Bregman"20G2VM arkevitch et alJ2003(l . Some 
striking visualizations of the kinetic SZ fluc tuations in- 
duced by these internal flo ws are given in Nag ai et al.l 
II2003H It has been shown iHaehnelt fc Tegmarklll996l: 
iHoldeil 120031 iNagai et alJ I2003D that the average pecu- 
liar velocity provides an unbiased estimate of the true 
bulk velocity but with an added dispersion of roughly 
100 kms~^. For our purposes this can be considered as 
an extra source of noise which is small compared to the 
uncertainties due to astrophysical confusion. 

We have neglected calibration uncertainty. Since the 
total signal, anywhere other than near the thermal SZ 
null, is about ten times as large as the kinetic SZ signal, 
accurate calibration of one band relative to another is 
important. Calibration errors will have to be controlled 
to better than about 10% for peculiar velocity measure- 
ments to be better than about one a. The same argument 
applies to gas temperatures for clusters with tempera- 
tures near 6 keV, since the relativistic correction is also 
about 10% of the total signal. But gas temperatures, 
since they do not suffer from CMB anisotropy confusion, 
can be measured to much better than 1 a. Achieving 
±1 keV, as can be done if we neglect calibration uncer- 
tainty, may require calibration uncertainties as small as 
2%. Measurements closer to the thermal SZ null will 
reduce these sensitivities to calibration error. 

6. EXPERIMENTS 

As our reference experiment we have modeled the sen- 
sitivity achievable for a generic ground-based bolometric 
experiment with three frequency bands selected to lie in 
the atmospheric windows available at a good millimeter 
site. We have also assumed that a 30 GHz measurement 
is available, with sensitivity similar to that expected from 
the SZA. We assume that the sensitivity in each channel 
of the bolometric experiment is limited only by fluctu- 
ations in the photon background itself and that other 
potential sources of noise, such as phonon fluctuations in 
the bolometers or electronic readout noise, are negligible 
(this can be achieved with current detector technology, 
see Lange 2002). We assume that 1// noise introduced 
by fluctuations in atmospheric water vapor emission can 
be adequately subtracted from all of the data. This can 
be accomplished by some kind of spatial chopping or re- 
moval of common-mode signals from array data, or alter- 
natively the spectral properties of the atmospheric noise 
can be used. The latter technique has been demonstrated 
by the Sunyaev-Zeldovich Infrared Experiment (SuZIE) 
which operates at 150, 220 and 350 GHz and which has 
been used to set limits to the peculiar velocities of 11 



galaxv clusters ijBenson et al.l2003albfl . For the purposes 
of this paper we do not include in our calculations any of 
the overhead (increase in real observation time) that this 
process will produce, since this will depend on the indi- 
vidual experiments. For example, since the atmosphere 
is in the near field of most large telescopes, many detec- 
tors view the same column of atmosphere, yet observe 
different parts of the cluster, or none of the cluster at 
all. By using detectors that are a large distance from the 
cluster to subtract common-mode fluctuations, it may 
be possible to subtract atmospheric noise with little or 
no degradation in sensitivity. This is a commonly-made 
assumption that we also include here. 

With the above assumptions the performance at each 
frequency is characterized by the noise equivalent power 
(NEP; quoted in WHz^^/^), given by: 

NEP2 = 2P^,^Ahv + (32) 

where Pioad is the total background loading on the detec- 
tor in Watts, Ai/ is the detection bandwidth in GHz and 
n is the number of waveguide modes that are detected 
(n — \ for a diffraction-limited system). The first term 
in equation 1321 is caused by shot noise due to Poisson 
statistics of the incoming photons while the second term 
accounts for the effect of photon (boson) correlation. For 
more detail, see iLamarrel ()1986j) . Sources of background 
loading include the atmosphere, warm emission from the 
telescope and surroundings, and of course the CMB it- 
self. The total background loading is then: 

Pload =Y.^^ (33) 

i 

where the sum is over all sources of power on the detector. 
For each source of loading, the power at the detector can 
be calculated from: 

P.= P-^^^^^mA^,Av (34) 

where Trj^; is the equivalent Rayleigh- Jeans (RJ) tem- 
perature of the ith background source, ?/; is the fraction 
of photons from the background source that reach the 
detector, A is the effective telescope area and fib is the 
solid angle response of each detector on the sky. The 
integral is over the spectral band to which the detector 
responds. We assume for simplicity, a square band of 
width Az/. We now consider each source of background 
power in turn. 
For the telescope the equivalent RJ temperature Trj.i 

is: 

Trj,! = e(i/)Ttoi— (35) 

where Ttei is the physical temperature of the telescope, e 
is a frequency-dependent emissivity and xj^ = hv /kT'^^y. 
For the atmosphere: 

Trj,2 = Tatm {1 - exp[-r(z/)/cos^za]} (36) 

where Tatm is the physical temperature of the telescope, 
riy) is the frequency-dependent zenith optical depth and 
6'za is the zenith angle of the observation. 
For the CMB: 

Tk3,3 = Tcmb— — - (37) 
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where Tcmb — l.ll^Y^ and x — hv jkTcMB- 
The optical efficiency is assumed to be: 
_ / 0.4 for i = 1, 2 (Tel, Atm) 

~ j 0.4 X exp[-T(i/)/ cos6'za] for i = 3 (CMB) 

which are values typical of those measured for bolomctric 
systems (Holzaofel et al. 1997; Runvan ct al. 2003). 

Once the NEP is determined, the sensitivity to CMB 
fluctuations, including the SZ effect, is obtained by calcu- 
lating the noise equivalent temperature (NET) as follows: 

NET - ^ K s-i (39) 

'\/2?7cMB ("^^cmb/^Tcmb) 

where ?7cmb is given above. 

We further assume that the focal plane is fully sam- 
pled with detectors spaced at half the beam size. Each 
detector will have an angular response that is diffraction- 
limited at each frequency, in which case the solid angle 
response of each detector is just fib = /{Av"^). With a 
suitable scanning strategy, an oversampled map can al- 
ways be produced with a pixel solid angle fJp. The rms 
sensitivity per map pixel is then: 



NET 



(40) 



where iint is the integration time per detector and the 
NET is per focal plane element. The factor / accounts 
for the increase in sensitivity that a fully sampled fo- 
cal plane can achieve over a focal plane with feed horns 
that are exactly diffraction-limited. We assume / = 1.58 
which is the ideal limit. The real value of / depends 
upon factors suc h as detector noise and instrument pho- 
ton background l)GrifBn et al.ll2002|) . 

We include a map pixel size here to allow us to com- 
pare experiments with different angular resolutions. If 
the array is large enough the scanning can be arranged 
so that there are always focal plane elements viewing the 
cluster. No integration time is lost in this case. Other 
chopping schemes could degrade the sensitivity by a fac- 
tor that depends on the details of the particular scheme 
being used. Thus our final numbers may be taken as best 
case values. 

In order to use eauationsl32lthrough l40l we now define 
the details of the experiments. For our reference exper- 
iment we assume a ground-based off-axis 8-m telescope 
located at a site with an altitude of 17,000' and precip- 
itable water vapor of 0.5 mm. The telescope and instru- 
ment are assumed to have an equivalent Rayleigh- Jeans 
temperature of Trj^i = lOK that is frequency indepen- 
dent. Figure 121 shows the dependence of 5T on frequency, 
for a map pixel size of 1 sq arcmin. Two models are 
shown, one with a bandwidth of 20% at each frequency 
and one with a bandwidth of 10%. There are broad at- 
mospheric windows visible at 135-165 GHz and at 215- 
290 GHz. The window at 345 GHz is sufficiently narrow 
that only a 10% wide frequency band can be utilized 
but it is a very useful atmosphere and/or point source 
monitor. Based on this figure we have placed three of 
the four channels of the reference experiment at 150, 220 
and 280 GHz. The sensitivities are shown in Tabled 

The Table also shows sensitivities for several experi- 
ments cur rently being built . The Sunyaev-Zeldovich Ar- 
ray fSZA. iCarlstrom et alJ 1)2002 )) is a 30 GHz interfer- 
ometer which will begin operation in early 2004. The 
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Fig. 2. — Error in each arcminute pixel for an hour of integration 
with our reference experiment as a function of central frequency, 
for bandwidths of 20% (solid line) and 10% (dashed line). 



SZA will also be able to observe clusters with higher an- 
gular resolution, but lower sensitivity, at 90 GHz (not 
shown in the Table). SuZIE HI will operate at the Cal- 
tech Submillimetcr Observatory (CSO) and will observe 
a 6' X 2' field of view simultaneously at 150, 220, 280 and 
350 GHz with four 48-pixel arrays. SuZIE HI has some- 
what worse sensitivity than the reference experiment be- 
cause the CSO is an on-axis Cassegrain telescope, leading 
to higher telescope loading due to the secondary mirror 
support legs obstructing the beam. The Atacama Cos- 
mology Telescope (ACT) is a proposed 6m telescope to 
be sited on the Atacama Plateau in Chile. The ACT will 
map one hundred square degrees of sky simultaneously 
in three frequency bands with 32 x 32 element arrays and 
is expected to detect large numbers of clusters through 
the SZ effect. 

7. RESULTS 

Here we display and discuss our error forecasts. Many 
of our plots show error forecasts as a function of 
the galaxy cluster parameters and Tg. We use 
the core radius as a proxy for redshift; for a fixed 
three-dimensional gas distribution the central y value 
is redshift-independent, so the redshift-dependence of 
the signal comes entirely from oc 1/Da{z) where 
Da{z) is the angular diameter distance to the cluster. 
The core radius is also an interesting parameter since 
it has a large influence on the peculiar velocity error 
l|Haehnelt fc Tegmarkl lT996^. Cluster mass is related to 
gas temperature by M = 8xlO^'^h-'^MQ{T/6keV)^^'^ (al- 
though the normalization is not well-established) so the 
temperature range 3 to 15 keV corresponds to a mass 
range 3 x lO^-* to 3 x 10^^ /i-^Mq. 

We study first the dependence of our error forecasts 
on galaxy cluster parameters and point source parame- 
ters. Then we look at the dependence on experimental 
parameters. Finally, we forecast errors for some planned 
experiments. We vary many parameters and each time 
hold many other parameters fixed. Table El shows the 
values of parameters we use unless otherwise specified. 
For the reader's convenience, this is a redundant list; for 
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1.0 


10.0 


SuZIE-III'' 


150 


1.0 


9.2 




220 


0.7 


15 




280 


0.5 


28 




345'^ 


0.5 


110 


ACT« 


145 


1.7 


3.4 (3.4) 




225 


1.1 


3.7 (3.4) 




265 


0.9 


4.3 (3.4) 



Resolutions assumed for our forecasting, for simplicity, 
are 1' for all experiments and all channels, except for ACT 
where we assume 1.7' for all channels. 

^ The sensitivity numbers are per 1' pixel assuming / = 1.58. 

Sensitivity for mapping an individual cluster in 12 hours 
(Leitch private communication) 

Sensitivity achievable in one hour over a 48 sq. arcmin 
map. 

^ This channel is unused in our forecasting. We assume it is 

used entirely for atmospheric subtraction. 

^Numbers in parentheses are from 

http://www.hep.upenn.edu/angelica/act/act.html for a 

100 sq. degree map. These are the numbers we use in our 

forecasting. 

Table 1 Assumed experimental specifications and 
sensitivities per square arcminute pixel. 



example is specified as well as T even though O is 
simply fcr/(meC^). 

7.1. Effect of Point Sources 

We begin with our fiducial cluster, studying the error 
forecasts as a function of confusion noise and spectral in- 
dex rms in Fig. |21 For the lower (Tq values one can infer 
from the flattening of the error contours with increas- 
ing shot noise amplitude that the data themselves are 
capable of fitting for the dusty galaxy component with 
residuals at about the 100 to 200 fi K level. If instru- 
ment noise were to increase, this transition from depen- 
dence on (^(^TfgQ^ )^/^ to independence would occur at 
higher values. The trend of error level with CTq, is as ex- 
pected: fiat at low ([ST^^f^] )^/^ and steeper at higher 

values. Peculiar velocities are more affected, particularly 
at non-zero CTq, than the other parameters. 

We now fix the point source contamination parameters 
and study our results as a function of galaxy cluster core 
radius, 9c, and gas temperature, T. Note that as we vary 
9c the map size we assume increases in order to assure 
complete coverage of the cluster and good CMB sub- 
traction. For large-format bolometer arrays, expected to 
have more than 1000 elements, the larger map size does 
not require longer integration times because the field of 
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Fig. 3. — Forecasted uncertainties for one hour of observation 
with the Reference 1', 4- band experiment and fiducial cluster pa- 
rameter values of dc = 30" and T = 6keV as a function of confusion 
noise rms in units of arcminutes and spectral index rms in units of 
/iK-arcmin. Fiducial point source parameter values are aa = 0.3' 

and ((5T|?o) )V2 = i7o^K-arcmin. 



view will be much larger than the galaxy cluster. There- 
fore, despite the added coverage, one hour of observation 
is still sufficient to obtain the Reference experiment sen- 
sitivities in Table 1. On the left of Fig. 21 we set the 
point source contamination parameters to zero and on 
the right to their fiducial values CTq, = 0.3/.tK-arcmin and 
(((5T3f^)2)i/2 = l70^K-arcmin. 

Because of the linear response of the signal to y, a{y) 
is independent of y and therefore T as seen in the left 
side of Fig.^l In contrast, at fixed y and r/3 the signal is 
proportional to so one expects cr{T) cx l/{dT^/dT) oc 
1/T. Confusion with y, more important at high T, leads 
to the error decreasing more slowly than 1/T. The cr(T) 
has a similar dependence on T since the negligible error 
on y means <t(t)It = a(T)/T and our scaling assumes 

TCX.T. 

For fixed three-dimensional distribution of gas pres- 
sure, central y values do not vary as a function of redshift. 
However, core radii do as 9c oc 1/Da{z) where Da{z) 
is the angular-diameter distance to redshift z. Thus to 
study dependence on z we vary 9c with fixed y. Note that 
Da{z)/Da{0.5) ~ 0.5, 0.7, 1.3, 1.4 for z = 0.2, 0.3, 1, 
2 respectively, assuming = 0.3 and JIa = 0.7. Note 
that the number of pixels with SZ signal in any given 
range scales as 9^, so the errors on r, y and T decrease 
as l/9c- 

The peculiar velocities on the other hand suffer from 
contamination from the CMB. As 9c increases it becomes 
more difficult to distinguish the signal proportiona] to t(3 
from the very red CMB power spectrum, hence a{f3) goes 
up. At large 9c, where the error in /3 is dominated by the 
CMB-induced error in t/?, the error in r is negligible. In 
this case, a{(3) = (1/t)(t(t/9) so a{(3) oc 1/T, which is the 
scaling we see. At smaller 9c uncertainty in r is no longer 
negligible and the scaling with T is more complicated. 

Including the effect of the fiducial point source con- 
tamination (right side of Fig. Q does not alter the error 
forecasts dramatically. For v the largest impact is at the 
smaller 9c — which is unfortunate since this is where v 
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Fig. 4. — Forecasted uncertainties for our Reference 1', 4-band experiment as a function of gas temperature and core radius. On the 
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is measured best. At small 6c the CMB noise is smaller 
and (less importantly) the Icnsing enhancement of con- 
fusion noise is largest. At 6c — 30" the velocity errors 
are almost doubled. 

7.2. Effect of Varying Experimental Parameters 
Figure [S] shows how the error forecasts depend on ob- 
serving time. On the left side one can see <j{y) oc yJ\/tohs 
indicating that for the Reference experiment the domi- 
nant source of uncertainty on y is instrument noise, not 
dusty galaxies. Similar scalings are seen for cr(T) and 
cr(T). The a{v) in contrast, quickly saturate with in- 
creased observing time producing very little improve- 
ment since the dominant source of error is the CMB 
contamination. The saturation occurs even sooner for 
larger 6c. 

As can be seen in Fig. parameters other than v are 
highly independent of the angular resolution Oh- This 
is because none of the signal is lost with increasing 0f,: 
the signal is simply spread out over greater area. Errors 
would go up with increasing 6b for single-pixel observa- 
tions as opposed to the map-making observations we con- 
sider. The peculiar velocity dependence on Oh is steep for 
^f) ^ 26'c, since spreading out the signal degrades the abil- 
ity to subtract off the CMB contamination. Note that to 
explore the effect of angular resolution independent from 
sensitivity, we have not scaled the sensitivity with beam 
size as one would expect from Eq. 1401 but rather fix the 
weight per solid angle. The apparent decrease in a{v) 
with increasing 9b at fixed 9c ~ 45" is an artifact of our 
sparse sampling of the 9c, 6b space. 

We now consider the importance of each of the chan- 
nels of our reference experiment. On the left side of 
Fig. Owe show error forecasts for the Reference experi- 
ment with no 30 GHz channel, 'Reference-30'. Without 
this channel it is harder to distinguish a change in y from 



a change in T, thus the errors in both of these increase 
considerably. Their ratio, r = y / {T /m), suffers similarly. 
The impact on peculiar velocities is less dramatic, espe- 
cially at high 9c where the CMB contamination remains 
dominant. Changes are significant at smaller 9c where, 
e.g., a{v) goes from 250 to 400 km/sec at T = 6 keV. 

The strong degeneracy between y and T for Reference- 
30 can be seen on the right side of Fig. [3 The correla- 
tion coefficient, rxy = {6TSy) / {a{T)a{y)) increases from 
0.93 for Reference to 0.998 for Reference-30. This near- 
maximal correlation coefficient means that the error in 
T (y) is increased by a factor of 14 compared to the error 
in T (y) if y (T) were held fixed. Including a reasonable 
prior on the velocity of 300 km/sec decreases the T and y 
errors by about 20%. The importance of a low frequency 
measurement for determining T has been emphasized by 
Holder ( 2003) andpighanim et al. ( 2003). 

After the 30 GHz channel, the most critical channels 
in decreasing order are 150 GHz, 220 and 280. Removal 
of the 150 and 220 GHz channels can increase a-{T) by 
up to a factor of 5, although only about a factor of 2 
at our small fiducial core radius of 6c = 30" assumed in 
the right side of Fig. [7\ These increases, at the assumed 
level of instrument noise, are not sufficient to degrade 
the peculiar velocity errors. 

The dashed Hue on the right side of Fig.[7|is y (x T^, 
through our fiducial values of y and T. We see that the 
degeneracy direction is such that assuming the scaling 
relation and a normalization would lead to much more 
precise determinations of both y and T. On the other 
hand, the orientation of the degeneracy makes testing 
the scaling relation more difficult than if it lay along the 
line. 



7.3. Planned Experiments 
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We saw that for the Reference experiment the loss of 
the 30 GHz channel increased cr(T), as well as a{v) for 
small Oc- The lack of a low frequency channel can there- 
fore be compensated by an X-ray determination of the 
temperature. Thus we show forecasts for the ACT sur- 
vey and SuZIE-III observations both complemented by 
X-ray temperature determinations of ±1 keV. The re- 
sults for the two different experiments are very similar. 
We assume one hour of integration with SuZIE-III for ev- 
ery 48 sq. arcmin covered. Recall that in order to assure 
good subtraction of the CMB, we assume a map of linear 
extent 9m — 4(20c + Ob)- We have made no attempt to 
optimize this observing strategy. 

The right-most figure in Figure |S1 shows that adding 
a 30 GHz measurement to the millimeter data can break 
this degeneracy, as we expect from our previous forecasts 
for the Reference experiment in Fig. 01 In 1 hour of ob- 
servation with SuZIE III (per 48 sq. arcmin of map) 
and a 12-hour SZA observation, the peculiar velocity of 
a cluster can be measured to 250 km s~^ even without 
an X-ray measurement of temperature. A low frequency 
channel will thus be crucial for obtaining peculiar veloc- 



ity and gas temperatures at high redshift where X-ray 
brightness is very low. 

We also see in Figure |H1 a temperature error for our 
fiducial cluster of ~ 1.5 keV (and even smaller for larger 
9c)- At this accuracy level, one expects to begin to see 
discrepancies between X-ray temperature measurements 
which are more sensitive to the hotter, inner regions and 
SZ temperature measureme nts which are more sensitiv e 
to the cooler, outer regions ijMathiesen fc Evrardll200lfl . 
Deeper integrations can slowly improve these measure- 
ments as seen from the study of dependence of forecasts 
for the Reference experiment on observation time. 

8. RADIO POINT SOURCE CONTAMINATION 

We have assumed that (i) the frequency spectrum of 
cluster radio sources will be sufRciently steep at high 
frequencies that measurements at 150 GHz and higher 
will not be contaminated and (ii) that measurements at 
30 GHz will be sufficiently high resolution that point 
sources can be subtracted at that frequency (this is pos- 
sible with interferometric experiments such as the SZA) . 
Consequently we have concentrated only on the contam- 
ination introduced by dusty galaxies known to be strong 
submillimeter- wavelength emitters and have excluded ra- 
dio sources from our analysis. We will now justify this 
omission. 

Clusters often contain radio point sources at mJy 
levels, as seen from 30 GHz SZ measurements 
(|LaRoaue et al.ll2003t iReese et all 1200^ . It is impor- 
tant to note that radio point sources and submillime- 
ter point sources are not the same objects. The ex- 
cess of radio sources toward galaxy clusters is most 
likely due to emissi on from cluster galaxies themselves 
ijCoorav et alJIlQQSl) , whereas the dusty galaxies are typ- 
ically background sources that are not associated with 
cluster members. Very little is known about the emis- 
sion from radio point sources at millimeter wavelengths, 
but their spectra are expected to steepen due to the 
energy losses of the most energetic electrons to syn- 
chotron emission. This is confirmed to some extent by 
spectral indices measured b etween 10 -90 GHz for bright 
sources ijHerbig fc Readhead,1992.: Sokasian et al. 2001) . 
although there is substantial scatter in the measured in- 



11 



Reference-50 a(T) = 10 keV 




10 10 2 4 6 8 10 12 14 

T/keV T/keV 
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GHz channel (not shown) makes the least difference of all. Dashed line in right panel shows expected scaling of McCarthy et al 2003. 
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Fig. 8. — Forecasted uncertainties for planned experiments as a function of gas temperature and core radius. Left and center are for 
ACT and a one-hour SuZIE-III observation respectively each with a ±1 keV determination of the temperature from X-ray data. Right is 
for a combined twelve-hour SZA observation with a one-hour per 48 sq. arcmin SuZIE-HI observation. 



dex. iSokasian et alJ 1)200 ll) find a mean index of —0.5 be- 
tween 10 and 90 GHz, while Hc rbig fc ReadheadI l|1992D 

find a mean index at 40 GHz of — 0.8±0.4. Measurements 
by Trushkin ( 2003) of the spectral indices of t he 208 point 
sourc es detected by WMAP at 40 GHz (Be nnett et al.l 
120031) . suggest a mean of [3 — —0.1 at 20 GHz with large 
scatter, but also show that the spectra of radio sources 
are often not well-approximated by a simple power law. 

Sources of only a few mJy at 30 GHz could easily still 
have mJy fluxes at 150 GHz if the spectrum does not 
steepen significantly above 90 GHz. We have used the 
point source c ounts from WMAP, D ASI, VSA and CBI 
summarized in lBennett et al.l l)2003j) to determine a best 
fit model of: 



dN 
dSv 



Sq 



where Nq = 30 per sq. deg. and 5*0 = 1 mJy. Our expres- 
sion gives number counts for mJy sourc es that are slightly 
lower than lWhite fc MaiumdaiH2003D. but that a r e con- 
sistent with those determined by lLaR,oaue etal\ l|2002fl 
based on SZ measurements of clusters at 30 GHz. Conse- 
quently we expect that the number of sources above 0.1 
mJy at 30 GHz is less than one per 4 sq arcmin (roughly 
the size of a moderate cluster). Repeating our Fisher 
analysis with radio sources at this level, we find that for 
all the experimental situations we have considered, the 



peculiar velocity uncertainty for our canonical source is 
unaffected, even if a spectral index of /3 = 0.0 is assumed, 
and the number counts are 10 times higher than we have 
assumed. 

The spectra of mJy radio sources will soon be much 
better understood when the SZA and the Array for Mi- 
crowave Background Anisotropy (AMiBA) start opera- 
tion at 90 GHz. However, a complete understanding of 
faint radio sources at mm wavelengths will likely require 
the sensitivity and resolution of the Atacama Large Mil- 
hmeter Array (ALMA) ^ 

9. APPLICATION OF GAS TEMPERATURE AND PECULIAR 
VELOCITY MEASUREMENTS 

At roughly the level of 1 keV the distinc- 
tion between X-ray emission-weighted temperatures 
and electron-weighted temp eratures becomes important 
(jMathiesen fc Evrardll200"H) . so a comparison of the two 
"gas temperatures" could be valuable for studies of the 
physics of the intra-cluster medium. 

Furthermore, the pressure-weighted temperature 
should be an excellent indicator of cluster mass. This 
could allow independent calibration of the mass scale 
for galaxy cluster surveys that use X-ray, optical or SZ 
selection criteria. 

^ http://www.alma.nrao.edu 
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The measurement of pecuhar velocities using the ki- 
netic SZ effect has been a goal of SZ measurements for 
many years. Peculiar velocities at large distances will 
allow measurements of cosmological parameters and re- 
construction of the large scale gravitational potential. 

The two-point function of radial peculiar velocities is 
sensitive mainly to Qm and largely insensitive to the 
dark energy equation-of-state parameter, Wx or other 
cosmological parameters. A linear theory calculation of 
how w ell the two-point function can constrain gives 
(IPeeTfc Knox .2002) 



0.04V4007/V 



(500 km s-i)" 



(42) 



The SuZlE 111 and ACT collaborations plan to measure 
optical redshifts for 400 clusters with masses greater than 
3 X lO^^A/0. Since typical errors for peculiar velocities on 
these may be 300 km/sec and the velocity rms is about 
400 km/sec, these data may be able to achieve Aflm — 
0.04. 

However, Equation 1421 assumes that the variance in 
our peculiar velocity measurement errors is perfectly well 
known. In reality, the statistical properties of the errors 
will be difficult to know well because of their dependence 
on the subtraction of the contamination of dusty galaxies. 
If the assumed error variance is different from the actual 
error variance by /(w^) then there will be a systematic 
error in fJ,,, of 



9 In (w^ 



0.02(//0.1) (43) 



independent of the number of clusters measured. To con- 
trol this systematic error, we will need to learn more 
about the spectral properties and luminosity functions of 
sub-mm galaxies. To avoid this source of systematic er- 
ror one could use a linear statistic, the average difference 
between the radial velo cities of two clusters separated by 
distance r ijJu szkicwicz et aDll999D . How this statistic 
depends on cosmological parameters is under investiga- 
tion. 

Another apphcation of pecuh ar velocities is gr avita- 
tional potential r econstruction l|Dekel et alJll990() . As 
iDore et all l|2003(l point out, due to the low number den- 
sity of clusters this is only possible on very large scales. 
Comparison with reconstruction from galaxy number 
counts can determine galaxy biasing properties. A good 
understanding of the large scale potential would be valu- 
able for investigating environmental effects on galaxy for- 



mation and would be a starting point for constrained re- 
alizations of simulations of the large-scale structure of 
our universe. 

10. DISCUSSION AND CONCLUSIONS 

Contamination by emission from dusty galaxies will 
make accurate measurements of peculiar velocities diffi- 
cult. With four observing frequencies, few fiK sensitiv- 
ity, and arcminute resolution it will be possible to mea- 
sure peculiar velocities at the level of roughly 200 km/s 
for massive galaxy clusters. 

For comparison, previous estimates of uncertainties in 
peculiar velocities by Holder (2003) assumed sensi- 
tivity and perfect point source removal, which on Figure 
^ roughly correspond to aa — and an observing time 
of several tens of hours. In that case it was found that 
uncertainties well below 100 km/s were possible, as could 
be extrapolated from Figure^ Such exposure times may 
be prohibitively long and removing point source contri- 
butions to this level will likely require ALMA. ALMA 
will provide both a better understanding of the physical 
nature of dusty galaxies in general as well as measure- 
ments of the contaminating galaxy fluxes at exactly the 
frequencies of interest. 

The Reference experiment, as well as combinations of 
planned observations, can provide measurements of pe- 
culiar velocities at the level of 150-200 km/s, even in the 
presence of contamination by dusty galaxies. This is suf- 
ficient accuracy to constrain to better than 10%. At 
the same time galaxy cluster temperature measurements 
will be possible at the keV level. This will allow inter- 
esting comparisons with X-ray spectroscopic measures of 
the X-ray emission-weighted temperature. 

Peculiar velocities measured using the kinetic SZ 
effect, a long-standing challenge to CMB experimental- 
ists, may soon be a reality. The final limiting obstacle 
appears to be contamination from dusty galaxies, but 
sensitive multi-frequency measurements should clear 
this hurdle. 
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